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A comprehensive data base of if - -atom level shifts and widths is re-analyzed in order to study 
the density dependence of the ^-nuclear optical potential. Significant departure from a t e frp form 
is found only for p(r) / po J$ 0.2 and extrapolation to nuclear-matter density po yields an attractive 
potential, about 170 MeV deep. Partial restoration of chiral symmetry compatible with pionic 
atoms and low-energy pion-nuclear data plays no role at the relevant low-density regime, but this 
effect is not ruled out at densities of order po and beyond. J?-nuclear bound states are generated 
across the periodic table self consistently, using a relativistic mean-field model Lagrangian which 
couples the K to the scalar and vector meson fields mediating the nuclear interactions. The reduced 
phase space available for K absorption from these bound states is taken into account by adding 
an energy-dependent imaginary term which underlies the corresponding J^-nuclear level widths, 
with a strength required by fits to the atomic data. Substantial polarization of the core nucleus 
is found for light nuclei, and the binding energies and widths calculated in this dynamical model 
differ appreciably from those calculated for a static nucleus. A wide range of binding energies is 
spanned by varying the K couplings to the meson fields. Our calculations provide a lower limit 
of Tjf = 50 ± 10 MeV on the width of nuclear bound states for K binding energy in the range 
Bjf ~ 100 — 200 MeV. Comments are made on the interpretation of the FINUDA experiment at 
DA$NE which claimed evidence for deeply bound K~pp states in light nuclei. 

PACS numbers: 13.75.Jz; 21.30.Fe; 25.80.Nv; 36.10.Gv 

Keywords: Kaonic atoms; K nuclear bound states; Density-dependent X-nucleus interaction; K nuclear 
relativistic mean-field calculations 




2 



I. INTRODUCTION 

The existence of a KN unstable bound state, the 1 = A(1405) 7r£ resonance, about 27 MeV below the K~p 
threshold was first discussed by Dalitz et al. |l| and shown to be compatible with strong coupled-channel vector- 
meson exchange s-wave interactions in the 1 = KN — ttT, system. This strongly attractive 1 = KN interaction, 

3, 



plus the moderately attractive 1=1 KN interaction 



|, suggest that the K -nuclear interaction is also strongly 



attractive. Indeed, the Born approximation for the -ftf-nucleus optical potential constructed from the leading-order 
Tomozawa- Weinberg (TW) vector coupling term of the chiral effective Lagrangian [jj, 

C = - ^ , (i) 

where / w ~ 93 MeV is the pseudoscalar meson decay constant, yields sizable attraction V^ t ~ —55 MeV for a 
central nuclear density p ~ po = 0.16 fm~ 3 . Iterating the TW term plus next-to-leading-order terms, within an in- 
medium coupled-channel approach constrained by the KN — 7r£ — 7rA data near the KN threshold, roughly doubles 
this A'-nucleus attraction to about —110 MeV at p . It is found in these calculations (e.g. |4j) that the A(1405) 
quickly dissolves in the nuclear medium at finite densities, well below po, so that the repulsive free-space scattering 
length clk- p becomes attractive, and together with the weakly density dependent attractive a K - n it yields a strongly 
attractive density-dependent Hp' A'-nucleus optical potential: 

V %(r) = -— &oWpM , (2) 

PKN 

where pk n is the KN reduced mass and the in- medium isoscalar KN scattering length bo (p) (positive for attraction 
in the sign convention common in this field) is defined by 



b o(p) = 7;( a K- P (p) + a-K-n(p)) , b (po) - 0.9 fm . (3) 

However, when Vj£ t is calculated self consistently, namely by including Vj£ t in the in-medium propagator used in the 
e qmti „„ _ V' m a ^ ,™ th e — K-,^s — „ „„ ly 
moderately attractive, with depth between 40 and 60 MeV [5j,|6j,|7j. We call such if-nucleus potentials 'shallow'. In 
contrast, strongly attractive ^T-nucleus potentials, of depth between 150 and 200 MeV, appear to provide the best fit 
to the comprehensive K~ atomic data which consist of level shifts and widths across the periodic table flflQ. We 
call such AT-nucleus potentials 'deep', noting that the fits provided by such 'deep' potentials to the K~ atomic data 
are by far superior to those provided by the 'shallow' potentials 7J. 

Another approach to the construction of Vj£ t in dense nuclear matter has been the relativistic-mean-field (RMF) 
approach where in addition to an attractive Lorentz- vector mean potential, similar in its origin t o Eq. Q, an attractive 
Lorentz-scalar mean potential associated with the poorly known sigma term is included 




This approach 
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too yields a strongly attractive V^ t , the precise extent of attraction depending on the phenomenological constraints 
imposed on the RMF calculation. Assuming dominance of the Lorentz- vector mean potential proportional to /~ 2 , 
one could test the implications of partial restoration of chiral symmetry, in the sense originally proposed for pionic 
atoms 3 , by replacing in Eq. / w by a density dependent decay constant /* according to 



.t: 2 = fl-^p , (4) 

mi 



where a is the pion-nucleon sigma term |l6( . This implies a specific density dependence of the optical potential : 

W = M0) (5) 
oyH) 1 - 0.046(7/3 v ' 

where the ttN a term is given in MeV and p is the nuclear density in fm -3 . This chirally motivated density dependence 



has proved successful in pionic atoms 



and p is trie nuclear density m tm . 1 his chirally motivated den 
QQQQQQ and in low energy pion scattering 3.Q 



23j,|24| where its effect 

can l>o separated from similar cll'ivis due U> i lw i hrrsljold rnorjsy depem leuce ol'tlic nixlorlying ;r.Y amplitudes 2ll 3- 
The A-nuclear interaction is also strongly absorptive, which arises from the available one-nucleon absorption reac- 
tions KN — > ttY with approximately 100 and 180 MeV energy release for the final hyperons (Y) E and A, respectively. 
Among the known hadron-nucleus interactions at low energies, only the TV-nucleus interaction is as strongly attractive 

n 

and absorptive as the A -nucleus interaction 25]. 

Recent experimental rep orts on candidates for A- nuclear deeply bound states in the range of binding energy 
B R ~ 100 - 200 MeV 3 27, 3 3 

again highlighted the quest on of how attractive the A-nucleus interaction is 
below the KN threshold. Obviously, the 'shallow'-type potentials cannot generate deeply bound nuclear states in 
the energy range ~ 100 — 200 MeV, whereas the 'deep' potentials might do. The distinction between 'deep' and 
'shallow' A-nucleus potentials becomes somewhat fuzzy within a dynamical approach which allows for polarization of 
the nucleus by the strong A-nucleus interaction. The depth of the A-nucleus potential becomes then state dependent, 
thus depending on the binding energy B^. Indeed, strong polarization of the nucleus by the A interaction was 
established in our recent RMF calculations for the light nuclei 12 C and 16 30]. Consequences of such polarization 
are a central topic of the present work. 

In Section 2 we re-analyze the existing K~ atomic data in order to gain information on the density dependence of 
the A-nuclear optical potential in the nuclear surface and to extrapolate smoothly to nuclear matter density. This 
procedure yields a 'deep' attractive potential of about 170 MeV at po- We also test whether starting from a 'shallow'- 
type potential [^, ^, 7 , and superimposing on it a density dependence of the form Eq. could improve the fit to 
the A _ -atom data as it appears to do in pionic atoms and in low-energy pion-nuclear scattering 

24) . For K~ atoms the answer is in the negative, perhaps because the nonperturbative aspect of the A(1405) 



dominance at the low density regime encountered in kaonic atoms invalidates the chiral perturbation counting rules 



[31 . This does not rule out a greater role for partial chiral restoration of a similar sort at higher densities, and for 
deeply bound K-nuclcar states, or for kaon condensation as discussed first by Brown et al. |32| . 

The case for A-nucleus bound states with ~ 100 MeV, generated by a 'deep' potential, was made by Akaishi 



and Yamazaki 



who fitted the underlying t matrix to a partial set of the available low-energy KN data. For the 



most recent publication that elaborates on their calculations for very light nuclear cores, see Rcf. 



341. In these few- 



body calculations, the strongly attractive A-nucleus interaction caused a substantial polarization, or rearrangement 
of the nuclear core. In our dynamical RMF calculations 3^|, the K and the nucleons interact through the exchange of 
scalar (a) and vector (w) meson fields which are treated in the mean-field approximation. The presence of K induces 
additional source terms in the equations of motion for the meson fields to which the K couples. This affects the scalar 
and vector potentials which enter the Dirac equation for nucleons and leads to rearrangement of the nuclear core. It 
also leads to a deeper A"-nuclear potential with respect to static calculations which assume a rigid nuclear core (an 
assumption that is justified for K~ and p atoms 30j,|35j). By successively allowing the K to polarize the nucleons, 
and the polarized nucleus to enhance the AT-nuclear interaction, this dynamical calculation is made self consistent. 
Similar TV-nucleus RMF calculations have been reported by the Frankfurt group 

In the present paper we extend the results of Ref. | 30| in order to study the dynamical aspects of the A-nucleus 
coupling across the periodic table, and in particular to study the widths expected for deeply bound states over a wide 
range of binding energies. This is done by scanning on the coupling constants of the K - meson fields. In Section 3 we 
outline the A-nuclcus RMF methodology used in this work and discuss its extension to describe absorptive interactions 
by which the AT-nuclear bound states acquire a width. The calculation of the width is also done dynamically, since its 
evaluation requires the nonstatic nuclear density and the binding energy which determines the reduced phase space 
available to the absorption reactions, and both of these quantities keep changing from one step of the iterative solution 
of the equations of motion to the next. The imaginary part of the potential affects also nonnegligibly the calculated 
binding energy. This important aspect of the A-nucleus dynamics is missing in the very recent RMF calculation by 
Zhong et al. 38]. In Section 4 we show and discuss the results of calculations across the periodic table. Special 
attention is paid to 12 C and other light nuclear targets for which a recent FINUDA experiment at DA$NE, Frascati, 
claimed evidence for deeply bound K~pp states |29|. We find large dynamical effects for light nuclei such as 12 C and 
we discuss these in detail. Section 5 summarizes the present work with conclusions and outlook. 
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II. KAONIC ATOMS 



A. Motivation and background 



As explained in the Introduction, the motivation for re-analysis of a comprehensive set of kaonic atoms data is two 
fold. First is the question of 'deep' vs. 'shallow' real A"-nucleus potential, in light of recent possible experimental 
evidence for the existence of deeply bound kaonic states whose binding energies exceed the depth of the shallow type 
of potential obtained from fits to kaonic atoms data. Such shallow potentials are obtained when fits to the data are 
made with a simple '£p' approach, as described below, and also when self-consistent approaches are used (HO,!^. 
However, if the deep variety of potential is confirmed, then the dependence of the KN interaction on the nuclear 
density becomes of prime concern. This density dependence is the second point which motivated the re-analysis of 
kaonic atoms data although there have not been any new experimental results on strong interaction effects in kaonic 
atoms since the early 1990's. 

The method adopted here is rather similar to our earlier work |8|, |9J, llflj , namely, performing global fits to a large 
set of data which covers the whole of the periodic table, using an optical potential within the Klein-Gordon equation, 



[A - 2^{B + Vopt + V c ) + {V c + B) 2 ] V- = (h = c=l), (6) 

where V c denotes the static Coulomb potential for the K due to the finite charge distribution of the nucleus, including 
the first-order a{Za) vacuum-polarization potential, \x is the JsT-nucleus reduced mass and B = -Bjj + ir^/2, with Bg 
and T x standing for the binding energy and width, respectively. The interaction of K with the nucleus is described 
here in terms of an optical potential V ov t which in the simplest Hp' form is given by 

2pV opt (r) = -Ml + —^-^)bo(pn(r)+P P (r)) , (7) 

where p n and p p are the neutron and proton density distributions normalized to the number of neutrons N and 
number of protons Z, respectively, A = N + Z, and ton is the mass of the nucleon. In the impulse approximation l tp' 
approach the parameter bo is equal to the ^T-nucleon isoscalar scattering length under the sign convention common 
in this field, otherwise this parameter may be regarded as 'effective' and its value is obtained from fits to the data. 
For A >> 1, Eq. Q reduces to Eq. J5J discussed above in the context of nuclear matter. An isovector term is 
not included in Eq. (JJJ as earlier analyses showed |l(J| that such a term is not required by fits to kaonic atoms data. 
Compared to the previous analyses we have now used a modified set of parameters for the nuclear densities. The 
density distribution of the protons is usually considered known as it is obtained from the nuclear charge distribution 
by unfolding the finite size of the charge of the proton. The neutron distributions are, however, generally not known 



to sufficient accuracy. Experience with pionic atoms |20j and with antiprotonic atoms |33] showed that the feature of 
neutron density distributions which is most relevant in determining strong interaction effects in exotic atoms is the 
radial extent, as represented for example by r n , the neutron density rms radius. Other features such as the detailed 
shape of the distribution have only minor effect. This is observed also for kaonic atoms where the dependence of the 
quality of fits on the densities of the neutrons is less than what is observed for pionic and antiprotonic atoms. The 
neutron densities used in the present work were of the two-parameter Fermi shape and of the 'skin' variety, as detailed 
m Ref. |3. The values of r„ were those that yielded the best fit to pionic and antiprotonic atoms data. 

With the potential Eq. J7J and varying the complex parameter bo we obtain the same results as before [lo| when 
the data span a range of nuclei from 7 Li to 238 U. The value of % 2 is marginally larger than the corresponding value 
found in [lOj because of the slightly different nuclear densities used in the two calculations. The depth of the real 
potential for a typical medium-weight to heavy nucleus is about 80 MeV. Recall that the corresponding value for the 



shallow-type potentials is close to 55 MeV 
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B. Density dependence 

The possibility of departure of the antikaon-nucleus potential from the simple 'ip' form had been studied over a 



decade ago 




10J and significant improvement in the fit to the data had been achieved. The empirical parameter 



bo was made density dependent (DD) by replacing bo in Eq. J7J) by 

b b + Bo [^] a (8) 
Po 

where po = 0.16 fm~ 3 . In this way it was possible to respect the 'low density limit' (when a > 0) by keeping bo 
fixed at its free KN value and varying the parameters Bo and a. In view of the importance of the modification of 
the KN interaction in nuclear matter we have addressed this question again in the present work in a different way. 

We loosely define in coordinate space an 'internal' region and an 'external' region by using the multiplicative 
functions F(r) in the former and [1 — F(r)] in the latter, where F(r) is defined as 

F(r) = — !— - (9) 

with x = (r — R x )/a x . Then clearly F(r) — > 1 when r < R x — 3 a x , which is the internal region. Likewise 
[1 — F(r)] — > 1 in the external region. Thus R x forms an approximate border between the internal and the external 
regions, and if R x is close to the nuclear surface then the two regions will correspond to the high density and low 
density regions of nuclear matter, respectively. We have therefore replaced the fixed bo in the tp potential Eq. Q by 

b -» Bo F(r) + b [1 - F(r)] . (10) 
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TABLE I: Parameter values from global fits to kaonic atoms data, for the models of Eqs. ffl- llUH . Additional parameters for 
the fit using Eq. (1101 are R x o = 1-30 ± 0.05 fm, 5 X = 0.8 ± 0.3 fm, a x = 0.4 fin. Underlined parameter values were held fixed in 
the fit process. \ 2 values are for 65 data points, bo and Bo are in units of fm. 



potential 


x 2 


Re&o 


Imfco 


ReB 


ImBo 


a 


Eq. {Hp') 


129.4 


0.62±0.05 


0.92±0.05 








Eq. EJ ('DD') 


103.3 


-0.15 


0.62 


1.62±0.07 


-0.02±0.02 


0.24±0.03 


Eq. {TOJ ('F') 


84.4 


-0.15 





1.42±0.05 


0.66±0.05 





Here the parameter bo represents the large-r interaction, and it can be held fixed, e.g. at the free KN value. The 
parameter Bq represents the interaction inside the nucleus. As stated above, this division into regions of high and low 
densities is meaningful provided R x is close to the radius of the nucleus and a x is of the order of 0.5 fm. To enable 
global fits to be made over the whole of the periodic table we parameterized R x as R x = R x q A 1 / 3 + S x and varied 
in the least-squares fit the values of Bq, R x q and S x while gridding on values of a x . The parameter Re&o was held 
fixed at its free KN value but the results depend very little on its precise value. 

Table [I] shows results of \ 2 fits to kaonic atoms data over the whole of the periodic table using the three potentials 
described above. For the third option, Eq. (|10fl . the values of the geometrical parameters for a x = 0.4 fm are found 
to be R x q = 1.30 ± 0.05 fm, 8 X — 0.8 ± 0.3 fm, implying that the modification of the free KN interaction takes place 
at radii somewhat outside of the nuclear 'half-density' radius. Changing the value of a x between 0.2 and 0.5 fm made 
no difference. The ansatz Eq. I|10f> was applied to both real and imaginary parts of Bq, otherwise the fit deteriorates. 

Figure ^ shows the real part of the J^-nucleus potential for 58 Ni and 208 Pb for the tp model and for the two 
density-dependent potentials described above. The difference between the relatively shallow tp potential and the deep 
potentials is clearly seen. Figure [21 shows similar results for the imaginary part of the potential. Although the two 
density-dependent potentials have very different parameterizations the resulting potentials are quite similar. 

In Figure are plotted the F(r) functions for 58 Ni and 208 Pb, not as functions of r but rather as junctionals of 
the local nuclear density p, which is of prime interest in the present work. We note that for these two examples, 
representing medium-weight and heavy nuclei, the two curves are very close to each other. The two vertical dotted 
lines represent 10% and 50% of a nuclear matter density of 0.16 fm~ 3 . It is evident that the potential departs from 
the p(r) dependence for density lower than about 20% of its central value. 

Finally, we have tried to fit the K~ atoms data with the density dependence given by Eq. (JSJ) as a possible signature 
of partial restoration of chiral symmetry in dense nuclear medium [15J. However, no improvement on the shallow tp 
potential for K~ atoms was achieved when starting from the self-consistent amplitudes of Ramos and Oset 6>\ using 
for a values around 50 MeV 16j . That is in line with Fig. |21 where it is seen empirically that the density dependence 
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FIG. 3: The modifying factor F, Eq. @, as a functional of p, see text. The vertical dotted lines indicate 10% and 50% of 
nuclear matter density po = 0.16 fm -3 . 

is effective only at fairly low densities, considerably lower than the ones reached in pionic atoms (~ 0.5 po). 



III. K NUCLEAR RMF METHODOLOGY 

The theoretical framework adopted in this work for calculating K nuclear bound states is the relativistic mean field 
(RMF) model for a system of nucleons and one K meson |s(J. In this model, the interactions between hadrons are 
mediated by the exchange of scalar and vector boson fields which are treated in the mean-field approximation. In 
the calculations reported in Ref. 3(J and extended here, we employed the standard RMF Lagrangian Cn for the 
description of the nucleon sector, with the linear parameterization L-HS of Horowitz and Serot 39J, as well as with 
several nonlinear parameterizations, but mostly NL-SH due to Sharma et al. |40j. These parameterizations give quite 
different values for the nuclear compressibility and, therefore, are expected to provide different rearrangement of the 
nuclear core due to the presence of K. The (anti)kaon interaction with the nuclear medium is incorporated by adding 



to Cn the Lagrangian density Ck 



C k =V* KW i K -m 2 K KK - a<yK m K aKK . (11) 



The covariant derivative X> M = + ig^K^)^ describes the coupling of the (anti)kaon to the vector meson u>. The 
vector field u) is then associated with a conserved current. The coupling of the (anti)kaon to the isovector p meson 
is here excluded. This is a good approximation for N — Z nuclei which holds for the majority of nuclei considered 
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here, but not necessarily for the heavier 208 Pb nucleus. We note in this respect that optical-potential fits as described 
in Sect. [H]for kaonic atoms, including data for nuclei with excess neutrons up to 238 U, found no need to introduce 
isovector components. 

Whereas extending the nuclear Lagrangian Cm by the Lagrangian Ck does not affect the original form of the 
corresponding Dirac equation for nucleons, the presence of K induces additional source terms in the equations of 
motion for the meson fields a and loq to which the K couples: 

(-A + ml) a = -g<j N p s - g a Km K KK + {-g 2 a 2 - g 3 a 3 ) , (12) 

(-A + rn 2 )^o = +guNpv -2guK(uK + 9l>kUo)KK , (13) 
where ljk is the K energy in the nuclear medium: 



= \j m 2 K + g aK m K a + p 2 K - g^K^o , 

and ps and pv denote the nuclear scalar and vector densities, respectively. The additional source terms due to the 
K in Eqs. H12(l - (|13(l affect the scalar and vector fields (potentials) which enter the Dirac equation for nucleons. This 
leads to the rearrangement, or polarization of the nuclear core in the presence of K. 

The equation of motion for K is derived from the Lagrangian Ck using standard techniques. In order to preserve 
the connection to previous studies of kaonic atoms, the corresponding Klein Gordon (KG) equation of motion for the 
K is used in the form of Eq. © which is rewritten here 

[A-2p(B sp - +V opt + V c ) + (V c + B s - p f]K = (ft = c = 1) , (14) 

where the superscript s.p. in J5 s p - = B^' + iT K /2 stands for the single-particle K binding energy which is equal to 
the K separation energy only in the static calculation, as elaborated below. Moreover, in this particular form of the 
KG equation, it admits a straightforward extention, by allowing V op t to become complex, to consider dynamically the 
width r^- of K nuclear bound states. This is a new and crucial element in our calculations. The real part of the K 
optical potential V^pt in Eq. i|14|) is then given by 

^-^[is-a-S^-JlH, 

p 2 niK zrriK 



where S — g a KO~ and V — g^pc^o are the scalar and vector potentials due to the a and uj mean fields, respectively. 
We note that ReV^pt is explicitly state dependent through the (1 — [B s ^ 
multiplies the vector potential. This factor was disregarded in our Letter 

arose only implicitly through the dynamical density dependence of the mean-field potentials S and V 



s-p- j_ v c ]/mK) energy-dependent factor which 
where the state dependence of ReV opt 
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The binding (separation) energy of K, B^, in the combined K nuclear system is given by 

B R =B(p)-B( A Z) . (16) 

Here B( A Z) is the total binding energy of the nucleus Z. It is to be noted that the total binding energy B(AZ) of 
the K nucleus, besides including nucleon single-particle energies, mean-field energies and center-of-mass energy, also 
includes the K single-particle binding energy B s ^?' from Eq. (|14fl and the additional meson mean-field contributions 
from the source terms due to (anti)kaons in the KG Eqs. (|12fl - l|13fl for the a and to meson fields. Whereas Bg — B s ^' 
for the static calculation, B s ^' > Bg for the dynamical calculation, and this latter difference defines the nuclear 
rearrangement energy which is related to the polarization of the nuclear core by the K. These statements are 
demonstrated in Section llVl 

Since the traditional RMF approach does not address the imaginary part of the potential, ImV^pt was taken in 
a phenomenological tp form, where its depth was fitted to the K~ atomic data and the nuclear density p was 
calculated within the RMF model. We emphasize that p in the present calculations is no longer a static nuclear 
density. Here it is a dynamical entity affected by the K, which is embedded in the nuclear medium and interacts 
with the nucleons via boson fields. The resulting increased nuclear density leads to increased widths, particularly for 
deeply bound states. On the other hand, the phase space available for the decay products is reduced for deeply bound 
states, which will act to decrease the calculated width. Thus, suppression factors multiplying ImV^ p t were introduced 
from phase-space considerations, taking into account the binding energy of the kaon for the initial decaying state, 
and assuming two-body final-state kinematics for the decay products. Two absorption channels were considered. The 
dominant one for absorption at rest is due to pionic conversion modes on a single nucleon: 

KN^irY,, ttA (-80%) , (17) 

with thresholds about 100 MeV and 180 MeV, respectively, below the KN total mass. The corresponding density- 
independent suppression factor is given by 



- Mf V [Ml - (m. + m Y r][Ml - (m Y - m^] @(Ml ~ m ' ~ ^ ' (18) 
where M i = m,K + mjy, Mj = Af i — Bg. The second channel is due to non-pionic absorption modes on two 
nucleons: 

KNN — » YN (-20%) , (19) 

with thresholds about = 140 MeV lower than the single-nucleon threshold. This second channel represents in our 
model all the multi-nucleon absorption modes which are not resolved by experiment. The corresponding suppression 
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factor is given by 



f M$ 2 [Mf - (m N + m Y y] [M| - (my - m N )*] 

h = -773- \/ 7772 7 i T277772 7 72T 9 M 2 - my - m N ) , (20) 



where M 2 = to_k~ + 2tojv ; M? = M 02 — -B^. The branching ratios (quoted above in parentheses) are known from 
bubble-chamber experiments 4l| . Although multi-nucleon absorption modes are often modeled to have a power-law 
p a (a > 1) density dependence, a linear density dependence (a = 1) may also arise in multi-step reaction mechanisms 
|42j. Since our comprehensive AT~-atom fits |8j are satisfied with a ~ 1, we here assume a = 1 which means that 
/ 2 too is independent of density. We comment below on the effect of a possible density dependence of / 2 , reflecting 
perhaps a p 2 dependence of the non-pionic decay mode Eq. (|19fl at high densities. 

Since E final states dominate both the pionic and non-pionic channels j4l|. the hyperon Y was here taken as Y — E. 
Allowing A hyperons would foremost add conversion width to K states bound in the region ~ 100— 180 Me V. For 
the combined suppression factor we assumed a mixture of 80% mesonic decay and 20% nonmesonic decay |4l|. i.e. 

/ = 0.8 h + 0.2 h . (21) 

In the calculations below, a residual value of / = 0.02 was assumed when both /1 and ji vanish. 

The coupled system of equations for nucleons and for the electromagnetic vector field Aq, for the p meson mean 
field, and for the mean fields a and uq Eqs. (|12fl - l|13fl above, as well as the KG Eq. (|14|) for K~, were solved self 
consistently using an iterative procedure. Obviously, the requirement of self-consistency is crucial for the proper 
evaluation of the dynamical effects of the K on the nuclear core and vice versa. We note that self consistency is 
not imposed here on the final-state hadrons which only enter through their on-shell masses used in the phase-space 
suppression factors given above. For the main 7rE decay channel it is likely that the attraction provided by the pion 
within a dynamical calculation [6| is largely cancelled by the nuclear repulsion deduced phenomenologically for E 
hyperons 



IV. RESULTS AND DISCUSSION 
A. Bulk properties of K~ nuclear bound states 

In this subsection we show and discuss results of calculations of 'bulk' properties for K~ nuclear bound states, such 
as level widths, average nuclear densities and rms radii, and single particle energies. The calculations cover a wide 
range of binding energies, in order to establish correlations between some of these properties and to study effects due 
to the nuclear polarization. The empirical values and for the RMF K coupling constants, as found from a 
fit to kaonic atom data (third and fourth lines in Table I of Ref. ^^|), were used as a starting point for calculations. 
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A full dynamical calculation was then made for K~ nuclear states starting from the static tp imaginary potential 
obtained from the atomic fit with Im&o = 0.62 fm from Table [i] while entering dynamically in the iteration cycles 
the resulting nuclear density and the suppression factor / as defined by Ea. (|21[l . This dynamical calculation, for the 
light nuclei 12 C and 16 0, led to a substantial increase of the calculated binding energy. For example, whereas the 
static calculation gave B K - = 132 MeV for the Is state in 12 C, the dynamical calculation gave B K - = 172 MeV 
for this same state. In this work, in order to produce different values of binding energies, we focus on a particular 
way of varying the depth of the real if - -nucleus potential by scaling down successively g a K from its initial value 
and, once it reaches zero, scaling down g U K too from its initial value gjQ until the K~ Is state became unbound. 
Obviously, the good global fit to the atomic data is inevitably lost once the coupling constants are varied in order to 
scan over a wide range of binding energies. The reverse procedure of first scaling down g U K from its initial value g^ K 
to zero, and then reducing g a x from its initial value g^ until the K~ Is state became unbound, gave very similar 
results except for weak binding where the calculated widths are larger than 100 MeV and hence are outside of our 
direct interest. Furthermore, in order to scan the region of large values of Bk- , of order 200 MeV, we also scaled up 

we have shown that similar results are obtained 
for different starting values for g a x and g U K , adjusted to a 'shallow' variety potential [g. We have also shown there 
the little sensitivity to the nuclear RMF model version chosen, the linear version L-HS |39j and the nonlinear version 



g a K from its initial value g~^ K while keeping g U K = 9^k- In Ref. 



NL-SH |40j. Below we discuss some aspects of the sensitivity to the nonlinear version used, particularly for high 
nuclear densities. 

Figure0]shows calculated widths T K - as function of the binding energy B K - for Is states in r-- 2 C, ^^Ca, A: 208 Pb, 
for the nonlinear NL-SH version |40j of the RMF model. The dotted line shows the static 'nuclear-matter' limit 

. f r (0) 



rV" ; _ , (22) 



m K 

where / is the phase-space suppression factor of Eq. (|21Jl and Tj9- is given by 

r<?_ = — Im&o Po , (23) 

PKN 

for the static value Im&o = 0.62 fm used in the calculations and for po = 0.16 fm~ 3 . Eq. 123|) holds for a K~ 
Schrocdinger wavefunction which is completely localized within the nuclear central-density pa region. The additional 
factor (1 — B K - /ijik)^ 1 in Eq. (1221) follows from using the KG equation rather than the Schroedinger equation. A 
small correction term of order V c /mx was neglected in this factor. It is clearly seen that the dependence of the width 
of the K~ nuclear state on its binding energy follows the shape of the dotted line for the static nuclear-matter limit 
of r^- , Eq. (|22p. This dependence is due primarily to the binding-energy dependence of the suppression factor / 
which falls off rapidly until B K - ~ 100 MeV, where the dominant KN — + 7rE gets switched off, and then stays rather 
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FIG. 4: Dynamically calculated widths of the Is K -nuclear state in X - 2 C, x i°Ca and K ?° 8 Pb as function of the K binding 



y 



energy, for the NL-SH RMF model 40]. For K l Pb, results using the L-HS RMF model |39j| are also shown, by open squares 
The dotted line is for a static nuclear-matter calculation with po = 0.16 fm -3 . 

flat in the range B K - ~ 100 — 200 MeV where the width is dominated by the two-nucleon absorption modes Eq. 119fl . 
The larger values of width for the lighter nuclei are due to the dynamical nature of the RMF calculation, whereby 
the nuclear density is increased by the polarization effect of the K~ as shown in the next figures. We note that the 
widths calculated in the range B K - - 100 - 200 MeV assume values of about 50 - 60 MeV for K 12 G, 35 - 45 MeV 
for K 40 Ca, and about 25 - 30 MeV for K 2m Pb in the range B K - ~ 100 - 160 MeV, decreasing gradually with A 
to the 'nuclear matter' dotted-line value of about 25 MeV. However, as is clearly observed in the figure, for „?? 8 Pb 
beginning at B K - ~ 160 MeV, a strong sensitivity to the nonlinear version of the RMF calculation develops whereby, 
in the NL-SH version [4(| used here, the 208 Pb s.p. levels undergo significant crossings, the maximal nuclear density 
increases substantially, and as a result the nuclear rearrangement energy becomes significantly larger. This does not 
necessarily affect other bulk properties, such as the average nuclear density discussed below, that do not show a similar 
sensitivity. For this reason we also demonstrate in Fig. 01 by open squares the results of using the linear version L-HS 
[3^ of the RMF calculation for K 208 Pb above 100 MeV. The nonlinear and linear versions give practically the same 
results below 160 MeV, and we note assuredly that the behavior of the calculated width as function of the binding 
energy for the linear version follows closely the 'nuclear matter' dotted line without any interruption around 160 MeV. 
The compressibility of the nucleus for linear versions is about twice that produced by nonlinear versions which are 
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FIG. 5: Dynamically calculated widths of the Is K -nuclear state in K _ Pb as function of the K binding ener gy, for the 
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nonlinear RMF models NL-SH 1 1(1 (solid squares) and NL-TM1 |4q| (solid circles), and for the linear model L-HS |39[ (open 
squares). The dotted line is for a static nuclear-matter calculation with po = 0.16 frn -3 . 



more realistic in this respect, and hence it is more difficult to polarize the nucleus within the linear version. 

In order to explore further the sensitivity to the nonlinear version of the RMF calculation, we show in Fig. 0the 
calculated widths for „?? 8 Pb using in addition to the NL-SH version employed above 40] also the NL-TM1 version 

n , 

due to Sugahara and Toki |46J. The results of using the linear L-HS version 39] are also shown for B K - > 100 MeV, 
as well as the dotted line for the static muclear-matter calculation specified above. The two nonlinear versions agree 
with each other very closely below 160 MeV but depart from each other above, with each giving rise to enhanced 
values of the width. The width evaluated using the NL-TM1 version appears to behave considerably more regu larly 
than the NL-SH version . We note that similar calculations using the NL1 version, which was used in Ref. 
completely break down for the high nuclear densities encountered in this kinematical region. In contrast, the linear 
version l-hs y , with a considerably higher value of nuclear compressibility, does not give rise to any enhancement 
of the calculated width and it practically coincides with the results of the nuclear-matter static calculation. Another 
comment is that in the region below 100 MeV, the static nuclear-matter calculation produces higher values of width 
than any of the RMF versions, since it assumes a perfect overlap of the K~ wavefunction with the central-density 
region. 

To conclude the estimate of K~ widths we comment that replacing p by p 2 for the density dependence of the 
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FIG. 6: Calculated average nuclear density p for X - 2 C, x - 6 0, K i°Ca and x l 08 Pb as function of the Is K~ binding energy, for 
the same nonlinear RMF model as in Fig. [I] 

non-pionic decay modes Eq. (|19fl is estimated to increase the above values of the width by 10 — 15 MeV. This estimate 
follows, again, from the increase of nuclear density with B K - noticed above. Switching on the 7rA decay mode would 
increase further this estimate by 5 — 10 MeV in the range B K - ~ 100 — 180 MeV. We therefore expect that the 
estimate T K - = 50 ± 10 MeV in the range B K - ~ 100 — 200 MeV provides a reasonable lower bound on the width 
expected for light and medium-weight nuclei in any realistic calculation. 

The next four figures exhibit various nuclear properties, calculated dynamically, for nuclei containing a K~ nuclear 
Is state. Figure El shows the calculated average nuclear density p = 4 J p 2 dr as a function of B K - for the same 
K~ nuclear Is states as in Fig. 0|and for Is states in ^-0 which is studied below. The average nuclear density p 
increases substantially in the light K~ nuclei, for the binding-energy range shown here, to values about 50% higher 
than for these nuclei in the absence of the K~ meson. The increase of the central nuclear densities is even bigger, as 
demonstrated in the next two figures, but is confined to a small region of order 1 fm from the origin. These features 
point out to a significant polarization of the nuclear core by the Is K~ , particularly in light nuclei. 

Figure [7| shows calculated nuclear densities of j^-G for several values of the Is K~ binding energies. The purely 
nuclear density, in absence of the K~ meson, is given by the dashed curve. The maximal value of the nuclear density 
is increased by up to 75% in the range of binding energies spanned in the figure, and the enhancement is close to 
uniform over the central 1 fm, decreasing gradually to zero by r = 2 fm which already marks the nuclear surface. In 
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FIG. 7: Calculated nuclear density p of K - 2 C for several Is K~ nuclear states with specified B K - values, using the nonlinear 
RMF model NL-SH 401 as in Fig. H The dashed curve stands for the 12 C density in the absence of the K meson. 
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FIG. 9: Nuclear rms radius and neutron single-particle energies for K -0 as function of the Is K binding energy, for the linear 
RMF model L-HS Q (open circles) and the nonlinear RMF model NL-SH 40] (solid circles) used in Fig. ^| 

this fairly small nucleus, the density is enhanced over a substantial portion of the nucleus. This is different than in 
heavier nuclei, as shown in Fig. |HI for the medium- weight ^!°Ca. Here, the enhancement of the maximal density (at 
the center of the nucleus) reaches almost a factor of two, but it subsides almost completely by r = 2 fm which is still 
well within the nuclear volume (the nuclear radius of 40 Ca is about 3.5 fm). This localization of the nuclear-density 
enhancement is due to the localized Is K~ density. As a result, the average nuclear density of K _ Ca shown in Fig. 
El is only weakly enhanced as function of Bk- ■ In heavier nuclei such as 208 Pb, the localization of the Is K~ density 
is as pronounced and the resultant nuclear density is enhanced by more than a factor of two over a small portion of 
the nucleus, confined again to the region r < 2 fm. 
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The upper and lower panels of Fig. [5] show the calculated nuclear rms radius and the Is and lp neutron single- 

- , foi 



particle energies E n , respectively, for K l ~0 as a function of Bk- , for two versions of the RMF calculation, the linear 



model L-HS 



K- 

(open circles) and the nonlinear model NL-SH |4(J (solid circles). The hgure shows clearly that the 
polarization effect of the Is K~ bound state on the Is nuclear bound state is particularly strong. The differences 
between the linear and nonlinear models reflect the different nuclear compressibility and the somewhat different 
nuclear sizes obtained in the two models. It is interesting to note that the increase in the nuclear rms radius of „ _ O 
for large values of Br- is the result of the reduced binding energy of the lpi/2 state, due to the increased spin-orbit 
term. Note also that as B K - approaches zero we do not recover the values inherent in static calculations for the 
various nuclear entities in Fig. ©and Fig. because at B K - — both ReV^ t and ImV^ t are still far from assuming 
zero values. 

The polarization effects noted here in both Fig. and Fig. [5] are somewhat larger than in Fig. 2 and Fig. 3 of 
our Letter publication |30| . This is due to the energy dependent suppressive factor in front of the vector potential in 
Eq. (|15|) which was disregarded there and which requires stronger couplings of the K in order to reproduce the same 
value of B K - as before. However, the effect of this suppressive factor on the calculated widths is minor as long as the 
width is plotted as function of the binding energy, as done here in Fig. 0|and in Fig. 1 of Ref. 3(]], mainly because the 
dominant factor in the calculation of the width is the phase-space suppression factor / which is directly determined 
by B K - . A welcome by-product of the present revision is that the 'deep' potential fitted in section [D] to the K~ 
atomic data produces K~ nuclear bound states in light nuclei in the range of interest, B K - ~ 100 — 200 MeV, where 
several claims for bound states have recently been made 



26, 27, 28, 291. This is discussed in the next subsection. 



B. Selected examples 

The discussion in this subsection is related to the FINUDA experiment |29j at DA<I>NE, Frascati, which recently 
suggested evidence for a K~pp bound state produced in K~ absorption at rest on unresolved combination of 6 Li, 7 Li 
and C targets. The experiment looked for back-to-back Ap pairs, assigned to the strong decay K~pp — > Ap which 
is a special case of Eq. (11911 . In addition to observing Ap pairs which correspond kinematically to K~ pp es sentially at 
rest, a peak in the invariant mass of Ap pairs was found that corresponds in the interpretation of Ref. [29j to a K~pp 
cluster bound by B K - pp — 115 ± 6 ± 4 MeV with a decay width r = 67 ± 14 ± 3 MeV. We comment that the mere 
observation of back-to-back Ap pairs from the decay of a K~pp cluster does not mean that this cluster is self-bound in 
the nucleus, in as much as the observation of back-to-back np pairs in the weak hypernuclear decay of A hypernuclei 
does not mean that Ap clusters are self-bound in the nuclear medium. Since the ^-nucleus potential is at least as 
strong as the A-nucleus potential, the most natural expectation is for X-nuclear bound states similar to the well 
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TABLE II: Binding energy B K - calculated statically and dynamically, nuclear rearrangement energy -RjJuci ano - width Y K 



calculated dynamically, for the Is K - nuclear bound state in K -Li, K ^G, K - 6 0, K i°Ca and K 2 _ 8 Pb, using coupling-constant 



ratios a a = QcK / = and = ^a'/sIk = 0.95 



(i) 



nucleus 



B 3 ^ (MeV) 



B° y _" (MeV) 



RtZ (MeV) 



idyn 



(MeV) 



6 Li 
12 ( 

40 Ca 

208p b 



2 C 
; 



28.1 
71.6 
66.0 
88.7 
107.3 



35.6 
88.7 
72.6 
93.3 
108.5 



17.0 
10.6 
7.6 
4.2 
1.3 



229.1 
94.9 

100.4 
51.6 
27.3 



accepted spectra of A-nuclear bound states over the periodic table beginning with A = 3. The interpretation of the 
FINUDA events in terms of a K~pp bound state requires that its binding energy is independent of the nuclear target 
off which back-to-back Ap pairs are emitted. At present the data do not allow to extract a statistically significant 
signal from individual targets used in the experiment. Since the constitution of targets in the FINUDA experiment 
favors 12 C, we will assume for the sake of argument that this observed peak corresponds to a Is or a lp K~ bound 
state in 12 C, or perhaps in n C or n B corresponding to production by the (K~,n) or (K~ ,p) reactions respectively. 
The consequences of this working hypothesis for the other targets will then be considered. Subtracting B pp — 27.2 
MeV from the binding energy 115 MeV assigned in Ref. [2^ to the K~pp-c\uster binding, the relevant K~ binding 
energy is B K - ~ 88 MeV. Regarding the width assigned by the FINUDA experiment to the bound-state signal, it 
could partly arise from Fermi-motion broadening which in the Fermi gas model yields the following contribution: 



./<(—) >-< — > = = — =20.3 MeV ( PF = 270 MeV/c) . (24) 
2 y ttln rriN 5 V 7m» 

This will reduce (quadratically) the assigned width of about 67 MeV down to 53 MeV. A rather broad distribution of 
back-to-back Ap events may also arise without invoking a bound state of any sort, due to the combination of Fermi 
motion of the initial pp pair and the nuclear final-state interactions of the emitted Ap pair, as shown by Magas et al. 
[47 1 after submission of the present work for publication. 

In Table HP we show results of static and dynamical calculations for Is states in ^-Li, „_ C, K ^0, K 4 -G& and 
if 208 Pb, choosing the coupling constants g a K and g^K such that B K - ~ 88 MeV for the Is ground state of K ^?C in 
the dynamical calculation. These coupling constants are not close to those required to produce the K~ atoms best 
fit, and the K~ potential determined here is considerably weaker than the potential of type DD or F in Fig. ^ A 
comparison between the very different binding energies calculated for ^-Li and ^-G leads to the conclusion that the 
peak observed in the FINUDA experiment cannot result from both Li and C. The relatively weak binding calculated 
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for K -Li implies, in agreement with Fig. 0]for other nuclei, a huge value of over 200 MeV for the width. The width 
of the K -C Is state is smaller, T = 95 MeV, but is larger than the width of the observed peak. These conclusions 
do not change upon slightly increasing g u K from 0.95 to 0.985, so that Br- ~ 88 MeV holds for K X -C instead of for 
K 12 G. We then find B K - = 88.6(95.9) MeV and Y K - = 96.9(66.5) MeV for A , n(12) C, while K 6 _U is considerably less 
bound. 

We note in Table [H] the relatively strong variation of the calculated binding energies between ^L 2 C and j}~-0 which 
is due to the difference between the underlying nuclear density input for these nuclei in the static RMF calculation. 
Choosing only one of these nuclei for the mass-number A systematics, the results of the table may be summarized as 
follows: 

• The K~ binding energy increases monotonically as a function of A. However, without the Coulomb potential 
V c it is almost constant from 12 C on due to saturation. 

• The width decreases monotonically with A, which is a corollary of the T(B K -) dependence in Fig. 



Except for ^--Li where the binding energy is relatively small, the difference between the binding energies 



calculated dynamically and statically is substantial in light nuclei, decreasing monotonically with A, and may 
be neglected only for very heavy nuclei. 

• The dynamical effect is very strong for X _ 2 C where the increase in the dynamically calculated Bk- is 17 MeV 
with respect to the static calculation. This increase gets larger with B K - , as demonstrated in the previous 
subsection for a a = a u = 1 where it amounts to 40 MeV. 

• The nuclear rearrangement energy, defined here as = B s SPl — B K - in the dynamical calculation, decreases 
monotonically with A. For a given value of A, however (but not demonstrated here), Rj^ increases with B K - . 

Except for K -Li, all the K~ nuclei demonstrated in Table UTI have additional bound states for the specified coupling 
constants. Thus, both f^-C and x - 6 have bound lp states, with SJT = 16.6, 30.6 MeV, respectively, but with 
very large widths: T = 158, 123 MeV, respectively. We therefore have an unfavorable situation of overlapping Is 
and lp levels. 

In Table ITTT1 we compare between dynamical calculations for ^-iLi and ^i 2 C, requiring that B K - ~ 88 MeV now 
holds for the lp instead of the Is state in K l 3C This requirement is motivated by the recent calculations of Yamagata 
ct al. 48] according to which the lp production cross section in [K~ , N) reactions is higher than for the Is state. The 
requirement B K - ~ 88 MeV for the lp state in ^L 2 C is satisfied using kaon coupling constants only slightly different 



than those found optimal for the K atomic fit of Rcf. 



The width of this lp state, 89 MeV, is close to the value 



95 MeV listed in Table UTI for the Is state constrained to the same binding energy. The Is states calculated in this case 
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TABLE III: Binding energy B K - and width T K - calculated statically and dynamically for the K - nuclear Is and lp states 
in ^-Li and ^-C, using q ct = 1.1 and a u — 1.0. 

nucleus nl B%* (MeV) T% & ± (MeV) B d ^ (MeV) r^i 1 (MeV) 

6 Li Is 75.1 79.0 146.5 74.5 

lp 4.1 70.7 7.4 184.4 

12 C Is 137.6 34.5 181.3 55.7 

lp 67.5 97.8 88.6 89.3 



for the two nuclei are narrower than the respective lp states, simply because their widths are suppressed stronger by 
the phase-space suppression factor /. The increase in Br- for these deeply bound Is states provided by going from 
the static calculation to the dynamical one is very substantial, with exceptionally large value of 71 MeV due to the 
strong polarization that the 6 Li nucleus undergoes due to a deeply bound K~ . We note that for very shallow bound 
states, as demonstrated here for the #--Li lp state and as will be shown in Fig. 1101 below for the ^-G Is state, there 
is little difference between the static and the dynamical calculations for -Br--. In contrast, the width T^- is very 
strongly enhanced in K 6 _ Li owing to the increase of the nuclear density provided by the dynamical calculation, from 
71 MeV in the static calculation to 184 MeV in the dynamical calculation for the lp state. Although the width of the 
Is state hardly changes, it would have decreased by roughly factor of two due to the phase-space suppression factor 
/ if the increase of the nuclear density did not offset this decrease. 

A comparison between the statically calculated (open circles) and the dynamically calculated (solid circles) B K - 
and r^-- for the Is state in ^-C is shown in Fig. as a function of the coupling-constant strenghts. It is clear 
that for Br- > 25 MeV the dynamical calculation gives higher binding than the static calculation does, with the 
binding-energy gain increasing monotonically with B K - , and this effect becomes important for B K - > 50 MeV. The 
effect on T K - is more intricate: the dynamically calculated width is larger than the width calculated statically, except 
in the region around B K - ~ 100 MeV where the kinematical phase-space suppression factor / of Eq. I|21|l decreases 
rapidly. In this region, the dynamically calculated T K - corresponds to binding energy B K - which is associated with a 
considerably more suppressed value of / than in the static calculation, and this effect wins over the width gained by the 
increased nuclear density. We note that the width calculated dynamically for the Is nuclear state in x i 2 C does not fall 
below 50 MeV for the range of variation shown, whereas the corresponding limiting value of the statically calculated 
width is about 35 MeV. Another feature shown in Fig. ^Jconcerns the effect of the imaginary potential on the binding 
energy: the dynamically calculated binding energy B K - when ImI4 p t is switched off is shown by the dotted line. It is 
clear that the absorptive potential ImI4pt acts repulsively and its inclusion leads to less binding, particularly at low 
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FIG. 10: Is K binding energy and width in K 1 2 C calculated statically (open circles) and dynamically (solid circles) for the 



nonlinear RMF model NL-SH 40] as function of the loK and oK coupling strengths: a u is varied in the left panels as indicated, 
with a a — 0, and ot a is varied in the right panels as indicated, with — 1. The dotted line shows the calculated binding 
energy when the absorptive K~ potential is switched off in the dynamical calculation. 

binding energies. In fact, without ImVo P t we would always have > B s ^ & } . The repulsive effect of ImT4 p t gets 

weaker with B K - , along with the action of the suppression factor /, and beginning with B K - ~ 100 MeV it hardly 
matters for the calculation of Br- whether or not ImI4 p t is included. 
Finally, we comment on the outcome of making g^K density dependent, 



9u>k{p) 



(25) 



1 - 0.046ap ' 

similar to the density dependence on the rhs in Eq. (0), which according to Ref. 1 1 5] mi^lii siumbif <• pari ial n si m at ion 
of chiral symmetry in dense medium. The starting parameters for the static calculation corresponded to a 'shallow' 
^-nucleus potential of depths about 60 MeV for both real and imaginary parts. For this choice ^- 2 C and K l -0 
have just one bound state (Is) each, with binding energy B = 5.3, 22.8 MeV, respectively. The introduction of 
the density dependence Eq. I|25|l within a dynamical calculation increases the binding energy to about 70 — 80 MeV. 
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However, the width of these states is large, close to 100 MeV. We conclude that although this ansatz for density 
dependence did not improve the ht to kaonic atoms data where the need to introduce phenomenological density 
dependence appears at low densities, its role in the higher nuclear density regime |li| and for generating nuclear K 
bound states cannot be ruled out at present. 



V. SUMMARY 

In this work we re-analyzed the data on strong-interaction effects in K~ atoms across the periodic table in order to 
study the density dependence of Vj£ t and its extrapolation into nuclear-matter densities. A departure from the t c sp 
dependence was found at the 20% range of nuclear-matter density po- Partial restoration of chiral symmetry, such 
as discussed in low energy pion-nuclear physics, was found to be ineffective at this low-density regime, but its role at 
higher densities cannot be ruled out. A smooth extrapolation into the nuclear interior gave a deep K nuclear potential 
V^p t which would allow the existence of K nuclear bound states even for a static calculation. We then extended the 
purely nuclear RMF model by coupling the K dynamically to the nucleus. Negligible polarization effects were found 
for atomic states, which confirms the optical-potential phenomenology of kaonic atoms in general, and as a valid 
starting point for the present study in particular. Binding energies and widths of K nuclear bound states were then 
calculated across the periodic table with the aim of placing lower limits on the widths expected for binding energies in 
the range of 100 — 200 MeV. Substantial polarization of the nucleus was found in light nuclei such as 12 C and 16 O for 
deeply bound K nuclear states, resulting in locally increased nuclear densities up to about 2po. An almost universal 
shape was found for the dependence of the K width on its binding energy. The widths are primarily suppressed by 
the reduced phase-space for absorption of K , and are enhanced by the increased density of the polarized nucleus. 
The present results already provide useful guidance for the interpretation of the recent FINUDA experimental results 
[29! by presenting some calculated dependence on the atomic mass number A and placing a lower limit Y ^ ~ 50 ± 10 
MeV on K states bound in the range ~ 100 — 200 MeV. Our calculated widths for a K~ Is or lp bound state in 
12 C, in the kinematical region of the FINUDA signal, are about 90 MeV which is somewhat larger than the reported 
width. We noted that a lp bound state in 12 C, rather than Is bound state, is compatible with the 'deep' type 
of A"~-nucleus potential fitted to the atomic data and perhaps also with deeply bound K~ Is states in the range 
B K ~ 150 - 200 MeV for lighter nuclear targets such as 4 He 2(1 III I- We have argued that if the FINUDA observed 
signal corresponds to a K bound state in 12 C, then it is unlikely to arise in Li targets, thus emphasizing the need to 
get statistically significant data from separate nuclear targets. For lighter nuclear targets, where the RMF approach 
becomes unreliable but where nuclear polarization effects are found much larger using few-body calculational methods 
[3^. we anticipate larger widths than 50 MeV for K deeply bound states in the range Bk ~ 100 — 200 MeV. 
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